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Abstract 

We numerically study two conservative two-dimensional maps, namely the baker 
map (whose Lyapunov exponent is known to be positive), and a typical one (ex- 
hibiting a vanishing Lyapunov exponent) chosen from the generalized shift family 
of maps introduced by C. Moore [Phys Rev Lett 64, 2354 (1990)] in the context 

of undecidability. We calculated the time evolution of the entropy Sg = — glj^ ' 

{Si = Sbg = ~^^iPi^^Pi)i and exhibited the dramatic effect introduced by 
numerical precision. Indeed, in spite of being area-preserving maps, they present, 
well after the initially concentrated ensemble has spread virtually all over the phase 
space, unexpected pseudo-attractors (fixed-point like for the baker map, and more 
complex structures for the Moore map). These pseudo-attractors, and the apparent 
time (partial) reversibility they provoke, gradually disappear for increasingly large 
precision. In the case of the Moore map, they are related to zero Lebesgue-measure 
effects associated with the frontiers existing in the definition of the map. In addition 
to the above, and consistently with the results by V. Latora and M. Baranger [Phys. 
Rev. Lett. 82, 520 (1999)], we find that the rate of the far-from-equilibrium entropy 
production of baker map, numerically coincides with the standard Kolmogorov-Sinai 
entropy of this strongly chaotic system. 
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1 Introduction 



Dynamical systems wliicli present exponential sensitivity to the initial condi- 
tions are called chaotic {also called strongly chaotic). Indeed, almost all orbits 
are unpredictable for any finite precision calculation, even if the evolution was 
purely deterministic. But dynamical systems at the edge of chaos are also 
somehow unpredictable, they typically exhibit a power-law sensitivity to the 
initial conditions, and are called weakly chaotic. For both cases, the sensitivity 
to the initial conditions is typically given P^PfS] by 

«"-JSf-o«-^ + (^-«)vi""-" (1) 



where Ax()f:) is the time-dependent discrepancy between two initially close 
trajectories. The coefficient Xg is a generalized Lyapunov exponent, and q is 
an index associated with the entropy [1] 

Sg^ (qen■,S, = SBG^-T.Pi^^P^)^ (2) 

5 i=i 



where BG stands for Boltzmann-Gibbs. This entropy is at the basis of nonex- 
tensive statistical mechanics, which has received many applications for com- 
plex systems. Eq. ([1]) recovers, in the q = 1 limit, the usual exponential diver- 
gence ^ = e'^^* (Ai being, say for the simple one- dimensional case, the standard 
Lyapunov exponent). If g = 0, we obtain the simple result ^ oc t. In general, 
g 7^ 1 yields a power-law dependence. 

A meaningful statistical description is possible even when the maximal Lya- 
punov exponent vanishes. This fact has been illustrated in the Casati-Prosen 
triangle map [5j, a mixing and ergodic conservative system which presents the 
extreme case of linear instability [5J. This map satisfies that, in the infinitely 
fine graining limit (i.e., W oo), the g-entropy increases linearly with time 
only for the value of the entropic index q = [6j; its slope is expected to coin- 
cide with the g-generalized Kolmogorov-Sinai entropy rate Kg [T]. Furthermore, 
this value is expected to coincide with the g-generalized Lyapunov coefficient; 
in other terms, a Pesin-like equality [7] is expected to hold for generic q as 
well. 

Conservative dynamical systems leading to entropic indices g 7^ 1 are certainly 
interesting. The case that we have just mentioned, namely the Casati-Prosen 
map (for which g = 0), is one such example. In the present paper we study 
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various aspects along this direction, possibly with q different from both unity 
and zero. 



We are interested in studying the entropy production^ i.e., the rate of increase of 
the g-entropy, of two-dimensional maps and, if possible, to connect it with the 
Kolmogorov- Sinai entropy rate [5] Ki, which is a property defined for (strongly) 
chaotic dynamical systems. In this paper we present numerical results for the 
non-disipative baker map, and for one (from now on referred to as Moore 
map) of the shift-like dynamical systems proposed by Moore [9]. The latter is 
expected to present a power-law sensitivity to the initial conditions, i.e., g < 1 
(possibly g 7^ 0, in contrast with the Casati-Prosen map). 

To evaluate the g-entropy of the system, we first partition the phase space into 
W » 1 little equal cells, we then choose one of these cells, and put within 
>> 1 random initial conditions. As time t evolves, the N points spread over 
the phase space in such a way that, at each time t, we have a set of numbers 
{iVj(t)} (I])^iiVj(t) = A^, Vt), so that Ni{t) is the number of points inside 
the i-th cell. Then, for each value of t, we can consider a set of probabilities 
{pi{t) = Ni{t)/N} to find A^j points in the i-th cell. For achieving a numeri- 
cally meaningful definition of the probabilities pi, the condition A^ >> ly has 
to be fulfilled; we typically consider A^ = lOW. At t = 0, all probabilities but 
one are zero, hence 5*^(0), calculated through Eq. ([2]), vanishes Vg. In other 
words, before the system starts to evolve, we know all the information regard- 
ing the occupancy of the phase space. As t evolves, information is lost and 
Sq{t) starts to increase. In all cases, Sg will be bounded by its corresponding 
equiprobability value {Sqjmax = ^i-g'^ ^ whose g = 1 limit case yields InW. 
The g-entropy production is then defined as 

Kg = lim lim lim . (3) 



In practice, we take values of (A^, W, t) large enough so that our numerical 
result for Kg becomes independent from them. 

Summarizing, we numerically study direct snapshots of the occupancy of the 
space phase, and the time evolution of Sg(t). From these, we shall present that 

(i) Ki = Ki for the baker map; 

(ii) numerical precision plays a relevant role in the phase space occupancy and 
time evolution of Sg{t) in both baker and Moore maps; 

(iii) the time evolution of Sg partially (but not completely) reflects the time 
evolution of the phase space occupancy in both baker and Moore maps; 

(iv) g 7^ 1 for the Moore map. 
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2 The conservative baker map 




Fig. 1. The non-disipative baker map. 

The non-disipative baker map is an area-preserving chaotic map from the unit 
square into itself ^JLj. The procedure (see Fig. 1) is equivalent to the following 
transformation: 

'(2xi,W2) (0<Xi<l/2) 
{xt+i,yt+i) = { (4) 
i2xt-l,iyt + l)/2) il/2<xt<l), 



where Xt and yt are the coordinates of each point on step t, evolving to Xt+i and 
yt+i on step t + 1. Baker map is a mathematical model of a mixing system. 
It exhibits a strong sensitivity dependence on initial conditions, thanks to 
the stretching in the x-direction. Its two Lyapunov exponents (Ai = In 2, 
A2 = — In 2) are opposite as required by the preserving area condition. 

To study the q-entropy production we first asume that phase space is parti- 
tioned into a grid ofW = 100 x 100 equal cells, and randomly choose = 10^ 
random initial conditions inside one randomly chosen cell. As time evolves the 
points spread and, taking account the set of probabilities into each cell, we 
calculate the g-entropy as a function of time t for different entropic indices 
q. In Fig. 2, we show the results we obtained working with double preci- 
sion. Observe a first stage during which the g-entropies grow, and a second 
stage during which the g-entropies remain practically constant (at the value 
Sq = liigW = , with 5*1 = \nW) for a relatively long time. After that, 
during a third stage, the g-entropies decrease, eventually vanishing (the fact 
that the entropy vanishes means that all points are now within a single cell, 
which by no means necessarily coincides with the initially occupied cell). It 
looks like if they followed a decreasing law virtually identical to the tempo- 
ral inversion of the first-stage increasing law: the map presents some sort of 
time reversibility. If we analyze the first stage, we observe that the value of 
g = 1 corresponds to a finite g-entropy production per unit time. For g < 1 
the curve is convex, while for g > 1 it is concave. This is due to the fact that, 
in chaotic systems like baker map, i.e., with strong (exponential) sensitivity 
to the initial conditions, the usual Boltzmann-Gibbs entropy (g = 1) is the 
appropriate one. Fixing this value of the entropic index, we consider series of 
thinner partitions in phase space. The resulting Si entropy is shown in Fig. 2 
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Fig. 2. Time evolution of the g-entropy for the non-disipative baker map, using 16 
digit calculations. The top figure shows Sq(t) for entropic indices q = 0.80, 0.85, 
0.90, 0.95, 1.00, 1.05, 1.10, 1.15, 1.20 (from top to bottom) when W = 10^ and 
A*" = 10^. The bottom figure shows the results corresponding to g = 1, for typical 
values of PV" (Ty = 10^, 4 X lO'', 16 x 10^; TV = l^W). Notice that the bounding 
value for the g = 1 entropy corresponds, in all cases, to equiprobability, i.e., hiW . 

for different values of W . Observe that, increasing W , the first stage remains 
for a longer time, which assures that the rate of g-entropy production for this 
far-from-equilibrium evolution depends neither on the number of points of 
the initial ensemble that we used, nor on W . If we perform a least-squares fit 
on these data, we obtain a linear regression index R = 0.99999998. The slope 
of this linear rise is 0.6928±0.0002, which, with great precision, coincides with 
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In 2, the positive Lyapunov exponent. We conclude that the rate of 5*1 produc- 
tion coincides, even for a finite size of coarse graining cells, with the standard 
Kolmogorov- Sinai entropy rate ki, which, according to Pesin theorem, equals 
the positive Lyapunov exponent. In other words, we verify 

= i^i = Ai . (5) 



This result reinforces those of Latora and Baranger [12113], in spite of the crit- 
icism of [13]. Indeed, Latora and Baranger showed, for several chaotic conser- 
vative dynamical systems, that the entropy evolution in far-from-equilibrium 
processes includes a stage during which the entropy linearly increases with 
time, the slope being numerically equal to the standard Kolmogorov- Sinai en- 
tropy rate. All these results are consistent with the possible existence of a 
g-generalized theorem along the lines of the Pesin identity. 



16-digit calculations (W=10x10): 

t=0 t=10 t=43 t=49 



t=51 t=55 t=60 t=72 




20-digit calculations (W=10x10): 

t=0 t=10 t=43 t=49 



t=53 t=69 t=75 t=85 



CO 3 



BAKER MAP 

N=1DW 
initial condition: (0,0) 
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Fig. 3. Numerical study of the baker map, with controlled fixed precision. The 
sequences of the top figure exhibit the evolution in phase space with a fixed precision. 
The corresponding curves for Si{t) are shown in the bottom figure. The evolution 
of Si{t) corresponding to a higher fixed precision experiment (40 digits )is shown 
as well; the time reversal of the entropy is not observed before t = 100. 
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The decreasing behavior of the g-entropy during the third stage needs fur- 
ther explanations. The decrease of entropy seemed to reflect a time reversal 
in the dynamic evolution. To elucidate this point, we performed iterations 
with controlled precision. The sequence of images in Fig. 3 shows that this 
reversal effect is gradually lost when precision is increased. Indeed, a roundoff- 
induced non-ergodic behavior emerges along the coordinate corresponding to 
the direction which exhibits strong sensitivity to the initial conditions. As a 
consequence, the g-entropy starts to decrease. This non-ergodic behavior even- 
tually leads to a final collapse onto the fixed point (0,0), which behaves as 
a pseudo-atractor. Consistently, the entropy turns back to be zero. Similarly, 
one must be aware that, on the plateau (i.e., the second stage), the entropy 
remains practically constant, in spite of the fact that it is possible to detect, 
in the evolution in phase space, gradual changes of phase space occupancies 
(see the W = 100 illustration in Fig. 3). If we increase the precision, we ob- 
serve that the second stage is longer. In other words, the system presents an 
equal-probability occupation for a longer time. We conclude that the entropy 
decrease hopefully disappears at the infinite-precision limit. This observation 
has an important practical consequence, namely that special attention has to 
be systematically dedicated to possible precision effects. This is analog to the 
results of Longa et al. [15J, who revealed the importance of the precision in 
order to understand the phenomenon of coalescence of the trajectories in a 
chaotic system. 



3 The (conservative) Moore map 
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Fig. 4. The Moore map. Alphabetic symbols are written on the cells to show how 
local dynamics evolves 

The generalized shift family of maps proposed by Moore [9J are a class of 
dynamical systems that posses some sort of undecidability, as compared with 
other low-dimensional chaotic systems PITU] . We are interested in studing a 
paradigmatic one within the Moore family. It is equivalent to the piecewise 
linear map shown in Fig. 4. When this map is recurrently applied, the area 
in phase space is conserved, while the corresponding shape keeps changing 
in time, becoming increasingly complicated. This map appears to be ergodic. 
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Frontiers definition ii 

(t=0) (t=1) (t=5.10') 




Frontiers definition iii 



(t=0) (t=1) (t=5.10') 





Fig. 5. Time evolution of three Moore maps (denoted by I, II and III) which differ 
just in the definition of the mapping of the frontiers. Top figures: Snapshots of the 
evolution in phase space (t = 0, t = 1 and i = 5 x 10^), when starting with points 
exclusively on the frontiers. In the t = 5 x 10^ squares, we have also indicated typical 
trajectories. Bottom figures: Time evolution of 5o.3 for maps I and II, starting with 
a set of initial conditions within a small cell {W = 10^; N = 10^). 
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16-digit calculations {W=10x10): 

t=0 t=200 t=1500 t=2000 t= 10000 




100-digit calculations (W=10x10): 

t=0 t=200 t=1500 t=2000 t=10000 
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Fig. 6. Numerical study of the Moore map I. Top figures: Evolution of occupancy 
in phase space. Bottom figure: Evolution of 5*0. i- 



possibly exhibits a Lyapunov exponent Ai = 0, and, presumably, the diver- 
gence of close initial conditions follows a power-law behavior [9j. When we 
consider a partition of W equal cells and select random initial conditions 
inside one random cell, the points spread much slower than they do on the 
baker map. More precisely, they spread, through a slow relaxation process, 
all over the phase space, each orbit appearing to gradually fill up the entire 
square. 

When we study the dynamic evolution corresponding to this map we observe 
that, in spite of being an area-preserving two-dimensional map, i. e. a con- 
servative map, it presents a sort of unexpected pseudo-attractors. The three 
temporal sequences in Fig. 5, correspond to three maps which only differ in 
the definition of the mapping of the frontiers existing in the piecewise linear 
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100 200 300 400 500 



t 

Fig. 7. Time dependence of Sq averaged over the 10% quick-best spreading cells, on 
the far-from-equilibrium regime, for typical values of q. 



Moore map. Different pseudo-attractors emerge for different mappings of the 
frontiers. The bottom figures refiect how the pseudo-attractors affect the long- 
time evolution of the respective g-entropies. If we define a map such that the 
area is preserved but the frontiers are dissipated, no pseudo-attractor emerges. 
This fact shows that the pseudo-attractors are due to zero Lebesgue-measure 
effects related to the frontiers. Consequently, in high-precision numerical stud- 
ies, the roundoff-induced long-time non-ergodic behavior tends to disappear. 

On calculations, with say sixteen digits, we observe that Sq initially grows. 
Then, its increasing rate slows down so as to tend towards the constant value 
that characterizes equiprobability. When precision is gradually lost along time, 
and non-ergodicity behavior emerges, Sq starts decreasing. The influences of 
the precision and of the grid size are illustrated in Fig. 6. 

The study of the g-entropy growth on the far-from equilibrium regime reveals 
important fluctuations. Fluctuations are usual at chaos thresholds and they 
requiere efficient and careful averaging over the initial conditions [16]. We 
checked that roundoff effects are not important in this regime and that they 
cancel over random averages. The optimal g-entropy average process consists 
in averaging over the quick-best spreading cells, which should be taken as 
independent initial conditions [16]. The g-entropy averages over the 10% quick- 
best spreading cells presents, for all values of q, intrinsic oscillations (see Fig. 7) 
depending neither on averages nor on precision. We verify that g = 1 does not 
produce, even discounting the intrinsic oscillations, a linear increase of Sq{t). 
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In other words, the usual Boltzmann-Gibbs entropy is not the appropriate one 
for simply describing this system. We consistently expect g < 1. However, its 
precise value elluded us because the oscillations make the numerical analysis 
quite complex. The determination of q remains as a future goal. 



4 Main results 

Let us summarize our main results: 

(i) The coincidence between the Kolmogorov-Sinai entropy rate ki and the 
entropy production Ki is numerically shown on a far-from-equilibrium con- 
servative chaotic system (baker map). 

(ii) Roundoff-induced pseudo-reversibility and pseudo-attractors are found on 
two dimensional chaotic conservative maps. 

(iii) A statistical entropy may only partially reflect the changes of distributions 
in the phase space for any finite coarse graining and any finite precision. 

(iv) In order to have a finite entropy production for Sq we need a value of 
q which is definitively smaller than unity, i.e., the Boltzmann-Gibbs entropy 
does not appear as the most adequate tool. Deeper studies are needed in order 
to establish whether another value of q can solve this problem. 

As we see, the entropic issues related with time evolution for deterministic 
nonlinear dynamical systems that are conservative are not very different from 
those corresponding to the dissipative ones. 
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